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Finite systems such as atomic nuclei present at phase transition specific features differ- 
ent from those observed at the thermodynamic limit. Several characteristic signals were 
found in samples of events resulting from heavy ion collisions at and above the Fermi en- 
ergy. The concomitant observation of different signatures of a liquid-gas phase transition 
"q ■ in nuclei on a given sample strongly supports the occurrence of this transition. 

1. INTRODUCTION 

A liquid-gas type phase transition is theoretically expected for nuclei owing to the 
analogy between the nucleon-nucleon interaction and the van der Waals force. Experi- 
mentally, indications of phase transition have been reported many years ago by groups 
studying multifragmentation. In the last few years progresses developed in the two do- 
mains: statistical physics built the concepts associated to the definition and signatures 
of phase transition in finite systems, among which atomic nuclei; on the experimental 
side new and performant 47r multidetectors allowed to study carefully selected samples of 
events instead of working on singles data. Among those the INDRA array was used to 
detect charged products formed in a wide variety of nuclear collisions, which bring nuclei 
in a state where phase transition will occur. The quality of the exclusive measurements 
permitted to search for different signatures of the transition, as it seems obvious that the 
concomitant observation of several of them will strongly reinforce the hypothesis that a 
phase transition has occurred. 

2. EXPERIMENTAL EVENT SAMPLES 

The results presented in this contribution were obtained in experiments performed 
with the 4tt array INDRA operating at GANIL or at GSI and mostly concern systems, or 
sources, with about 200-230 nucleons: "fused sources" from Xe+Sn, Ni+Au reactions, Au 
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or Ta quasi-projectiles. In order to analyse data in the framework of statistical physics, 
one must select well defined sources. In central collisions this was done by two techniques: 
either a selection of compact sources by taking events whose main direction of emission 
deviates by more than 60° from the beam direction (Xe+Sn system) pQ; or for the asym- 
metrical Ni+Au system a discriminant analysis method [2|. In the case of "fused sources" 
the experimental distributions of excitation energy are narrow. 

The selection is more delicate when one wants to isolate quasi-projectiles, because of 
the presence of direct and neck emissions. For symmetric systems, a first order separation 
is to identify as quasi-projectile (quasi-target) the ensemble of products with velocities 
parallel to the beam direction larger (smaller) than the centre of mass velocity. This 
method was used for the Au+Au data discussed later. The advantage of quasi-projectiles 
is that the whole possible range of excitation energy (0 - available energy) can be explored. 

In all tolerance of about 10% on the source charge is imposed. 

Any of these data samples constitutes a statistical ensemble, in the sense of a collection 
of events with similar properties. To apply statistical physics analyses one must be care- 
ful about the correct framework (microcanonical, canonical ensembles, with or without 
dynamical constraints ... ). 

3. SIGNALS OF PHASE TRANSITION : THEORETICAL BASES 

While at the thermodynamic limit first order phase transitions are characterized by 
discontinuities or divergences of observables, smoother situations are expected in finite 
systems, which exhibit anomalous curvatures in thermodynamical potentials. Saying that 
the entropy, S(E), presents a convex curvature at first order phase transition is fully 
equivalent to stating that the energy distribution is bimodal, meaning that at a given 
temperature correspond two different values of the energy. Another consequence is that, 
in a microcanonical framework, the caloric curve, T(E), should present a backbending for 
the energy domain corresponding to the phase transition, and therefore the heat capacity, 
c=dE/dT, will have a negative branch. However caloric curves may differ depending on 
the thermodynamical path followed by the system in the temperature-pressure-energy 
space, and may not present any backbending, even at phase transition. A more robust 
signal was found in the fluctuations of the kinetic part of the energy, directly connected to 
the heat capacity: very large fluctuations sign a negative value of c. See |3l4j for reviews. 

In this picture, one then expects to observe at the same time a bimodality of the 
energy and a negative heat capacity as signals of a first order phase transition in nuclei. 
Moreover as the system is then located in the mechanically unstable spinodal region of 
the phase diagram, additional information can be found in searching for characterization 
of the dynamics of the phase transition, which in nuclei may proceed through nucleation 
or spinodal decomposition. Many theoretical predictions being in favour of the latter [S], 
a third signal of phase transition might be the evidence of spinodal decomposition. 

On another hand the first indications of phase transition in the nineties were related 
to critical behaviors, such as the observation of a power law mass distribution for light 
products of multifragmentation jH]. At the thermodynamic limit, scalings are general 
properties of matter near the critical point and thus typical of second order phase tran- 
sitions. Recent studies have however shown that scaling properties were observed on a 
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much larger range when the considered system is finite: it is sufficient that the correlation 
length be close to the system size to observe scalings. In a lattice gas framework, for in- 
stance, it was demonstrated that Fisher scaling occurs when the system lies on a "critical 
line" which extends well inside the coexistence region [7j. For nuclei scalings can therefore 
also be considered as indications of the occurrence of a first order phase transition. 

4. SIGNALS OF PHASE TRANSITION : EXPERIMENTAL EVIDENCES 

Most of the data mentioned in this section are published, we will not illustrate them 
here but rather send the reader to the quoted references. 

4.1. Dynamics of the phase transition: spinodal decomposition 

Theoretical studies as well as semi-classical simulations of collisions showed that, in 
the course of the reaction, the system is driven to the spinodal region of the phase dia- 
gram |8l9j . If low enough densities are reached, the system will never recover to normal 
conditions, but rather break-up into pieces. The rupture is initiated by local density fluc- 
tuations which are amplified in an homogeneous way by the mean field, creating regular 
high density regions which will finally turn into equal size fragments, surrounded by nu- 
cleons and very small fragments (the gas). The signature of spinodal decomposition in 
infinite nuclear matter is thus the observation in the exit channel of fragments of a given 
size, linked to the wave length of the most unstable mode. Here again the picture is dif- 
ferent for nuclei, because of their finiteness. Firstly in this case several modes have equal 
characteristic times, leading to fluctuations in the size and multiplicity of the fragments 
from event to event jHj. Beating of modes and coalescence during fragment separation 
may break the symmetry of each partition. And last, as fragments are born hot, the 
de-excitation stage also alters the partitions. As a result the favoured partitions resulting 
from spinodal decomposition of nuclei might only survive with a very small probability, 
making them difficult to reveal. 

The method proposed to evidence enhanced equal-size fragment partitions uses higher 
order charge correlations ^U], each event being characterized by the average and the 
standard deviation of its fragment (Z>5) charge distribution. The delicate point in any 
correlation function is the construction of uncorrelated events. The initial method did not 
account for charge conservation; the INDRA collaboration proposed two methods to rem- 
edy this drawback, which lead to different results. Applied to semi-classical simulations of 
collisions between Xe and Sn at 32A MeV, where it is known that spinodal decomposition 
does occur, the modified event mixing method did not show any print of favoured parti- 
tions ^1], while the intrinsic probability method (IPM) enlightens an excess of about 1% 
of events where fragments have equal sizes. On experimental data the results of the two 
methods also differ, the event mixing method always leading to less positive conclusions. 
With the IPM, evidence of spinodal decomposition was found, at the level of 0.3-1% of 
the events, in central collisions between Xe and Sn at 32, 39 and 45A MeV, the signal 
vanishing at 50 A MeV |12] . It appeared at 52 A MeV for central Ni+Au collisions, but 
not at the lower energy of 32 A MeV |T3]. The thermal energy put in the systems showing 
prints of spinodal decomposition was evaluated (with the help of a statistical model) to 
5-7 A MeV [H]. 
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4.2. Bimodality 

By definition, at first order phase transition, the distribution of an order parameter, 
or of any related variable, should be bimodal in a canonical framework. This was looked 
for in Au quasi-projectiles from Au+Au collisions between 60 and 100A MeV. To mimic 
a canonical, or a least a gaussian ensemble [T3] , the sorting for the quasi-projectile was 
performed as a function of the transverse energy of the light charged particles (lcp, Z< 2) 
emitted backward of the centre of mass, E® 2 . This is justified by the high efficiency 
of INDRA for lcp in the entire phase space. At first order E® 2 quantifies the dissipated 
energy. For a specific bin of E^ 2 a bimodal distribution is found for the charge asymmetry 
of the two largest fragments of each event, (Z max — Z max _i) / (Z max + Z max _i) , when plotted 
versus Z maxi see examples in fig. |2(a)| and |2(b)| two groups of events clearly stand out, 
large asymmetries correspond to a big fragment associated with one or several very small 
ones (residue, or liquid type), while asymmetries close to are associated with events with 
more and smaller fragments (gas type). Outside this specific region, charge asymmetries 
are close to 1 for low E^( 2 and turn to for high E^( 2 [Tol - 

In the transition E{ X2 region, the properties of the two groups of events were investi- 
gated. The small asymmetry events have larger multiplicities of both fragments and lcp, 
indicating that the excitation energy of these events is larger than that of the other group. 
This is confirmed by calorimetry measurements: the energy distributions associated with 
the two groups of events are distinct, although overlapping [15]. If one admits that a 
selection of E® 2 fixes the temperature, then two values of the energy are correlated to it, 
which is the signature of a bimodal behaviour. 

4.3. Negative heat capacities 

The microcanonical heat capacity is a thermodynamic variable which has the specific 
property to present a negative branch in the transition zone for finite systems. Experi- 
mental determinations of the heat capacity rely on the measurements of kinetic energy 
fluctuations at freeze-out [3]. The heat capacity is given by the relation c = c|/(cfc — (< 
A > xaD/T 2 ), where Ck = 5 < Ek/A > /ST is the kinetic heat capacity, A and T 
the mass and temperature of the system and the kinetic energy fluctuation. The heat 
capacity is negative when the fluctuation is larger than the reference kinetic heat capacity. 

The method suffers however several drawbacks, some from lacks in the measured quan- 
tities (number and energy of neutrons, fragment masses), and others intrinsic to the 
method (reconstruction of the freeze-out configuration from the measured products). The 
influence of the implied hypotheses was carefully tested on event samples obtained in a 
statistical model, SMM [T7j. Interestingly it was demonstrated that the method tends to 
decrease the fluctuations, meaning that large fluctuations, and thus negative heat capac- 
ities, are not an artefact created by the imperfections of the measurement. 

First observed on Au quasi-projectiles from 35A MeV Au+Au semi-peripheral colli- 
sions [IB] , negative heat capacities were also seen in INDRA data for the same reactions 
at 60-100A MeV incident energies [T^]. The results of both experiments are in good agree- 
ment. In central collisions the end of the negative branch and the recovering to positive 
values were seen in Xe+Sn [20] and Ni+Au ^3] collisions. The ranges of excitation ener- 
gies where negative values of c were measured are in good enough agreement in all these 
cases, if one reminds that in central collisions some collective energy has to be added to 
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the potential energy part. 
4.4. Scalings 

Two types of scalings were studied in the INDRA data. Fisher scaling was already found 
in many experiments, and critical exponents were extracted. A-scaling was proposed more 
recently for determining the scenario of a second order phase transition. 

4.4.1. Fisher scaling 



11(A) 




(a) Excitation functions for different nu- 
clei produced in 39A MeV Ta+Au semi- 
peripheral collisions. Symbols are exper- 
imental data, lines show a fit with the 
Fisher formula. 

Figure 1. Fisher scaling 
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(b) Scaled cross sections for nuclei from boron to 
phosphorus produced in semi-peripheral (Ta+X) or 
central (Xe+Sn, Ni+Au) collisions. Points for the 
different system, which all lie on a single curve, are 
shifted for a better view. 



Fisher droplet model was originally proposed to study the morphology of fluids in 
thermal equilibrium in terms of an ideal gas of clusters coexisting with a liquid fraction. 
This model was recently applied to multifragmentation data, by considering all fragments 
but the largest as the gas phase, Z max being assimilated to the liquid part. The yield of 
a fragment of mass A reads: diV/dv4 = 1](A) = q A~ T exp((AAfi — coeA a )/T). In this 
expression, r and cr are critical exponents, A/i is the difference between the chemical 
potentials of the two phases, Co is the surface energy coefficient; e = (T c — T)/T c describes 
the distance of the actual to the critical temperature. This kind of scaling was found in 
many multifragmentation data, in hadron-nucleus as well as in nucleus-nucleus collisions; 
the agreement between data and theory often holds over orders of magnitude, and the 
critical exponents which are deduced are in good agreement with those expected for 
the liquid-gas universality class. In INDRA data, Fisher scaling was applied to fragments 
obtained in single sources as well as in quasi-projectiles [2T] . An example of the excitation 
functions for light fragments emitted by tantalum quasi-projectiles is given in fig. |l(a)[ 
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the Fisher formula well fits the experimental points over large ranges of masses and 
temperatures. Note that here the temperature is replaced by the square-root of the 
excitation energy, measured by calorimetry, according to the Fermi gas relation E* = aT 2 . 
In fig. |l(b)| are shown the scaling properties of the data obtained for different fused systems 
(Xe+Sn, Ni+Au) or Ta quasi-projectiles from four different reactions; once scaled, all data 
collapse onto a single curve - they have been shifted along the abscissa on the figure for 
better visualisation. For quasi-projectiles, and for fused sources formed at 32A MeV, the 
scaling extends below and above the critical temperature. This unexpected behaviour 
comes from the finite size of the systems. As shown in for finite systems the "Fisher 
critical temperature" obtained is not the true critical one if the system is in the coexistence 
region, and it varies with the density of the system. It happens for the data shown in 
fig. |l(b)| that the critical excitation energies found are in all cases close to 4. 5 A MeV |21j . 
If one remembers that the systems depicted in the figure also have similar masses, and 
thus may be represented by a single phase diagram, the unicity of the critical energy might 
be an indication that the density at phase transition (when the system multifragments) 
is similar for all systems, be they fused systems formed in violent collisions or quasi- 
projectiles from more peripheral collisions. 

The values found for the topological exponent r (2.1-2.4) and the one related to the 
surface-volume dimensionality ratio o (0.66) are in agreement with previous data, close 
to the ones of the liquid-gas universality class. 



4.4.2. Universal scaling 

The theory of universal fluctuations of order parameters provides information on the 
relationship between the phase transition of a system and the formation of clusters, with- 
out requiring a precise knowledge of the thermodynamical state of this system [22] • All 
information is supposed to be contained in the multiplicity and size of the clusters. For 
systems which exhibit second-order critical behaviours, the critical order parameters can 
be identified through their A-scaling behaviour. 

In the INDRA multifragmentation data on central collisions between a large variety 
of heavy ions, it appeared that the fragment multiplicity did not present any scaling, at 
variance with the size of the heaviest fragment, Z max . The nature of the order parameter 
identifies multifragmentation as an aggregation scenario [23j- Symmetric systems with 
total masses between 73 and 394 all showed two scaling regimes, from A=l/2 at low 
incident energy to A=l at higher energies, indicating the passage from an ordered phase 
to a disordered phase; the transition energy decreases when the systems grow heavier, 
(see |23|24j for details). It is interesting to note that the scaling behaviour is robust; for 
central collisions, it is rather independent of the precise event selection: the choice of 
compact sources j2H] leads to Z max distributions which are slightly narrower and shifted 
towards lower values than those obtained with an impact parameter selector [23] , but the 
transition energy is barely displaced. 

A-scaling was also evidenced in quasi-projectiles, provided that a selection is made 
via the bimodality order parameter. Events with large asymmetries scale with A=l/2 
whereas those with more equal fragments scale with A=l, see fig. |2(d)| 

The A-scaling signature of a phase transition is up to now the only one which has been 
observed without ambiguity for systems of total mass lighter than 150. 
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(a) Transverse energy distribution of lcp mea- 
sured backward of the cm. velocity, normalised 
to the incident energy per nucleon. The bars 
show the region where bimodality occurs. 



(b) Charge of the largest fragment vs charge 
asymmetry of the two largest fragments, Z asym . 
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(c) Kinetic heat capacity (stars, dotted line) and 
kinetic energy fluctuations (other symbols, full 
line.) Lines are to guide the eye. 



(d) A-scaling. Events with large charge asymme- 
try lie on the right branch (A=l/2), those with 
small asymmetry are on the left branch (A=l). 



Figure 2. Quasi-projectiles from 80A MeV Au on Au collisions. Part (b) is obtained for 
events inside the bars in (a). See ^H] for details, (from |19j ) 



5. SIMULTANEOUS OBSERVATION OF SEVERAL SIGNALS 

As said above, the hypothesis that atomic nuclei undergo a first order phase transition 
is strongly strengthened if several signatures of the transition are found on a given sample 
of events. Two examples will be presented in this section. 
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5.1. Quasi-projectiles from Au+Au reactions 

The first case, shown in fig |21 is that of quasi-projectiles (QP) produced in Au+Au 
reactions at 80A MeV measured with INDRA at GSI. As explained in 14.21 the properties 
of quasi-projectiles were followed after sorting on E^ 12 , the distribution of which is shown 
in fig. |2(a)[ Bimodality neatly occurs in the delimited E^ region, as shown in fig. |2(b)[ 

The A-scaling was searched for in all E®2 bins, and appears when separating the two 
groups of events. Events associated with large asymmetries (or Z max ) belong to the 
A=l/2 family, while those with small asymmetries correspond to A=l: the transition is 
found in the energy region where bimodality occurs (diamonds and crosses in fig. |2(d)[ ). 
Finally the negative heat capacity was calculated in a microcanonical framework inside 
each E®2 bin. A further selection on event compactness was added for this study [T9"] . 
The excitation energy of quasi-projectiles is determined by calorimetry. In figure |2(c)| 
are shown the kinetic energy fluctuations (solid line), versus the QP energy. Different 
symbols correspond to different bins of E^J. Fluctuations overcome the kinetic heat 
capacity (dotted line) in a small excitation energy domain of the QP (3-6^4 MeV), signing 
a negative value of the heat capacity; the largest fluctuations are obtained for events 
inside the E^ bin where bimodality occurs (crosses in fig. |2(c)| ). Note that fluctuations 
are maximum for an excitation energy of ~4.5A MeV, which is also the Fisher critical 
energy for similar (Ta) QP, as indicated in the previous section. 

Finally, for the selected sample of Au quasi-projectiles excited at about 4-5/1 MeV two 
signals characteristic of a first order phase transition are found at the same time. The 
observation of a scaling law also supports the existence of such a transition, because the 
system studied is finite. 



5.2. Single-source events from Xe+Sn reactions 

The second example, displayed in fig El concerns fused sources formed in central xenon 
on tin collisions at 25, 32, 39, 45 and 50A MeV. The event samples are in all cases obtained 
with the same selection (compact event shapes). Fig |3(a)| shows the excitation function 
of the excess of events with equal-sized fragments, signature of a spinodal decomposi- 
tion. The trace of the phenomenon is maximum at 39A MeV. The end of the negative 
branch and the restoring to a positive value of the heat capacity are seen for the 32 and 
39 A MeV samples (fig. |3(b)| ). From the Fisher scaling one finds a critical excitation en- 
ergy of 4.5A MeV, explored in the reaction at 32A MeV only; one expects a change in the 
multifragmentation pattern around this energy, which matches the energy where the ki- 
netic energy fluctuations are maximum [21]. Finally the transition from order to disorder 
marked by the change in A scaling also occurs between 32 and 39 A MeV(fig. g(dj. 

In summary, in events from central Xe on Sn collisions we have evidenced four signals 
indicating that the phase transition region is explored for incident energies between 32 and 
39A MeV. The transition energies found in the different analyses do not fully coincide; 
calorimetry gives about 6A MeV, while a determination through a comparison with the 
SMM model indicates 5 A MeV, close to the Fisher critical energy of 4. 5 A MeV. All values 
are however in agreement within the present precision of about 1A MeV that we have on 
calorimetry measurements. 
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(a) Excitation function for the (b) Energy distribution and heat capacity for collsions at 
spinodal signal between 32 and 32 and 39 A MeV (from |17j). 
50AMeV (from [I2j). 
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(c) Fisher scaling at 32, 39, 45 and (d) A-scaling. Points from right to left were obtained for 
50A MeV(from [21]). collisions at 25, 32, 39, 45 and 50A MeV. (from [22]). 

Figure 3. Central collisions between Xe and Sn from 25 to 50 A MeV 

6. Summary and prospects 

The two above examples clearly indicate that one can characterize a first order phase 
transition in nuclei with masses around 200, independently of the violence of the collision 
in which they were formed, and of the amount of radial expansion energy (larger in 
central than in peripheral collisions). Work is in progress to look for signatures of phase 
transition on lighter systems. To be more quantitative and do metrology, refinements of 
the methods are in progress, and experiments with larger statistical samples are necessary. 
In the future more information will be obtained when new accelerators will furnish very 
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exotic beams, permitting to explore the influence of the N/Z degree of freedom on phase 
transition. Calculations predict for instance the shrinking of the spinodal region when 
the number of neutron becomes very large (2H|. Ultimately, reliable information should 
be gained on the symmetry term of the equation of state. 
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